Strong and Weak Chaos in Networks of Semiconductor Lasers 
with Time-delayed Couplings 



Sven Heiligenthal. 1 '[2] Thomas Jiingling, 1 Otti D'Huys, 1 Diana A. Arroyo-Almanza, 2- ^] 
Miguel C. Soriano, 2 Ingo Fischer, 2 Ido Kanter, 3 and Wolfgang Kinzel 1 

1 Institute of Theoretical Physics, University of Wurzburg, 97074 Wurzburg, Germany 
2 Instituto de Fisica Interdisciplmar y Sistemas Cornplejos, IFISC (UIB-CSIC), 
Campus Universitat de les Illes Balears, 07122 Palma de Mallorca, Spain 
3 Department of Physics, Bar-Ran University, 52900 Ramat-Gan, Israel 
(Dated: October 9, 2012) 

Nonlinear networks with time-delayed couplings may show strong and weak chaos, depending on 
the scaling of their Lyapunov exponent with the delay time. We study strong and weak chaos for 
semiconductor lasers, either with time-delayed self-feedback or for small networks. We examine the 
dependence on the pump current and consider the question whether strong and weak chaos can be 
identified from the shape of the intensity trace, the auto-correlations and the external cavity modes. 
The concept of the sub- Lyapunov exponent Ao is generalized to the case of two time-scale separated 
delays in the system. We give the first experimental evidence of strong and weak chaos in a network 
of lasers which supports the sequence 'weak to strong to weak chaos' upon monotonically increasing 
the coupling strength. Finally, we discuss strong and weak chaos for networks with several distinct 
sub-Lyapunov exponents and comment on the dependence of the sub-Lyapunov exponent on the 
number of a laser's inputs in a network. 
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I. INTRODUCTION 

The cooperative behavior of a system of interacting 
units is of fundamental interest in nonlinear dynamics. 
Such complex networks have a wide range of interdis- 
ciplinary applications ranging from neural networks to 
coupled lasers [1-5 . Typically, these units interact by 
transmitting information about their state to their part- 
ners, and in many systems the transmission time is larger 
than the time scales of the individual units. For this rea- 
son, networks with time-delayed couplings are a focus of 
active research [5J [7] . 

Time-delayed feedback can produce dynamical insta- 
bilities which may lead to deterministic chaos [5MlO|. 
Even a scalar nonlinear differential equation with time- 
delayed feedback has an infinite-dimensional phase space 
which favors chaotic solutions. As an example, a sin- 
gle semiconductor laser often exhibits chaotic emission 
dynamics when its laser beam is reflected back into its 
cavity by an external mirror. Networks of nonlinear units 
may, similarly, become chaotic due to time-delayed cou- 
pling of the nodes 

For networks of identical units, chaos synchronization 
is possible. Even if the delay times are very long, the 
units may synchronize onto a common chaotic trajec- 
tory without time shift (zero-lag synchronization) |12f - 
RH] . Other kinds of synchronization are possible as well, 
like phase, achronal, anticipated and generalized synchro- 
nization. Chaos synchronization is being discussed in the 
context of secure communication |15l 116] . 
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Recently, two different kinds of chaos have been identi- 
fied for chaotic networks of time-continuous systems with 
time-delayed couplings: strong and weak chaos |17j . In 
the limit of large delay times, the maximal Lyapunov ex- 
ponent (LE) of the network saturates at a nonzero value 
for strong chaos, whereas it scales with the inverse delay 
time for weak chaos. A similar phenomenon has been 
reported in [9_ for time-discrete maps with delay. Only 
networks exhibiting weak chaos can synchronize to a com- 
mon chaotic trajectory. 

In this paper we extend these investigations on strong 
and weak chaos focusing on the dynamics of semicon- 
ductor lasers. In Sec. II we define the model and the 
linearized equations which determine the Lyapunov ex- 
ponents. In Sec. Ill a single unit with time-delayed feed- 
back is considered. An artificial sub-LE is defined which 
determines, whether the unit is in the strong or weak 
chaos phase. The scaling argument of [T7] is extended 
to derive the behavior of the maximal LE in the limit of 
large delay time. This system is investigated in Sec. IV 
for semiconductor lasers. Numerical simulations of the 
Lang-Kobayashi equations yield the transition from weak 
to strong chaos and back to weak chaos upon monoton- 
ically increasing the coupling strength. The scaling just 
at the transition is derived. Auto-correlations, spatial 
representations of the chaotic intensity and external cav- 
ity modes are calculated to investigate whether one can 
deduce the type of chaos from a single trajectory. In 
Sec. V networks of semiconductor lasers are considered. 
The stability of the synchronization manifold is related to 
the eigenvalue gap of the coupling matrix and the master 
stability function. The first experiment on semiconduc- 
tor lasers which supports the sequence 'weak to strong to 
weak chaos' with increasing coupling strength is demon- 
strated. Finally, networks with several distinct sub-LEs 
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and certain network patterns are investigated. 

II. THE LANG-KOBAYASHI EQUATIONS AND 
THEIR LYAPUNOV EXPONENTS 

We study the dynamics of coupled semiconductor 
lasers by means of numerical simulations, complemented 
by experiments. To this end we use the Lang-Kobayashi 
(LK) model, which compresses the complex physical pro- 
cesses of a laser in rate equations for only three classical 
degrees of freedom, namely the complex electric field £(t) 
and the normalized deviation of the charge carrier density 
n{t) with respect to the density at the lasing threshold. 
The rate equations for a single laser without self-feedback 
read 

£(t)= 1 ^G N n(t)S(t) 

n(t) = (p - 1) J th - 7 n(t) - [T + G N n(t)] \£ {t)\ 2 

A table of the parameters involved can be found in 
the appendix. Note that the phase of the electric field 
4> e (t) — arg£(t) in the model is defined to describe only 
the deviation from ujq t, which contains the high oscilla- 
tion frequency of the actual laser mode. Moreover, we 
have omitted intensity dependent nonlinear gain reduc- 
tion, which would be taken into account via an additional 
factor (1 + e\£\ 2 )- 1 with e ~ 1CT 7 for G N . All results 
presented in this paper, which are obtained with the sim- 
plified Eq. 0, coincide qualitatively with the results of 
the complete rate equations. The LEs are calculated by 
linearizing Eq. Q. The linearization reads 

S£(t) =^T^ g n [n(t) 6£(t) + £{t) Sn(t)} 

Sn(t)=-( 1 + G N \£{t)\ 2 )Sn(t) ^ 
- 2 [T + G N n(t)} ft {£* (t) 6£{t)} . 

We describe the system's state by the state vector x = 
(xi, X2, xs) T with x\ = X2 = and x% = n. 

Consequently, Eq. (JTJ reads x = F(x) and the lineariza- 
tion accordingly 8x = DF(x) <5x. In this notation we 
introduce time-delayed feedback on a single laser by 

x = F(x)+ ( xH(x r ), (3) 

where H(x) = cxp(— i cuo t)(x\, X2, 0) t describes linear 
coupling by reinjection of laser light into the cavity. The 
carrier density is not directly affected by the coupling. 
The coupling rate a denotes the feedback strength and 
depends on the reflectivity of the distant mirror as well 
as absorbers in the light path. The exponential term is 
a remainder of the transformation from optical phase to 
relative phase. The linearization of the laser equation 
with delayed feedback reads 

Sx = DF(x) 5x + a DH(x T ) 6x T (4) 

and can be used to calculate LEs for a given trajetory 
x(t), which is obtained by integration of Eq. |3|. The 



maximal LE A m is defined by 

A m =hmi Infirm (5) 

The norm ||a|| in our case is \f a\+ <z|, meaning that we 
use only the electric field components of the Lyapunov 
trajectory Sx(t). The choice of the norm does not in- 
fluence the final result, but may have strong influence 
on the required length of the trajectory to calculate an 
exponent with sufficiently high accuracy. 

Since the system under investigation is a delay system, 
we formally obtain infinitely many LEs depending on the 
initial conditions of the linear system. In practice, for 
random initial conditions after some transient time, the 
system relaxes onto the most unstable mode revealing 
the maximal exponent. In the following, we concentrate 
on the description and prediction of the dependence of 
A m on the feedback parameter a. 



III. STRONG AND WEAK CHAOS 

A. The sub-Lyapunov exponent and its 
experimental measurement 

Conditional LEs play an important role for synchro- 
nization in networks of coupled dynamical systems. More 
general, when a large system can be divided into two sub- 
systems forming a drive-response scheme (driver: T>, re- 
sponse: 1Z), it makes sense to define a conditional LE A-r 
of the driven system 7Z in order to classify the dynamics 
or to predict synchronization properties |18] . The con- 
ditional exponent is often referred to as 'sub-Lyapunov 
exponent' to stress the fact, that it characterizes a sub- 
system. If the sub-exponent is negative, 1Z is in a state 
of generalized synchronization with T>. Otherwise 7Z is 
autonomously chaotic, meaning that its state is not de- 
termined by the state of T>. 

The sign of a sub-LE can be measured experimentally. 
To this end a copy TZl of 1Z has to be created. 1Z' is 
then coupled to T> in the same way that 1Z is, meaning 
both units receive the same driving signal. If Ak < 0, 
then 1Z' and 1Z synchronize completely, otherwise not. 
This procedure is known as the Abarbanel test [TU] for 
generalized synchronization between T> and 1Z. 

A similar situation is present in case of a system with 
time-delayed feedback. Although quite against intuition, 
one can imagine partitioning of the delay system into 
non-delayed dynamical unit A = 7Z and transmission line 
A T = T>. The sub-Lyapunov exponent of the unit arises 
only from the instantaneous part of the equations of mo- 
tion, and has therefore been referred to as 'instantaneous 
Lyapunov exponent' |17j . However, in this paper we keep 
the notion 'sub-Lyapunov exponent' and denote it as Ao- 
Formally, from the equations of motion 

x = F(x) + crH(x r ) (6) 
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Figure 1. (a) Configuration to measure the difference be- 
tween strong and weak chaos, (b) Simulated isochronal cross- 
correlation C between two lasers coupled according to the 
setup presented in (a) for a delay time r = 10 ns and a pump 
current p = 1.02 vs coupling strength a. 



we obtain Ao by integrating the linearization 

5x = DF(x) <5x (7) 



An = lim - In 



t-s-oo t y ||<$Xo(£o 



l**b(*)ll 



(8) 



Alternatively, we can define Ao via the evolution operator 
of Eq. |7] 



5xo(*) = I/(t J *o)*x (*o), 



(9) 



with U(t ,t ) = 1. Then 
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A = max lim — \n(ei gl {U T (t, t ) U(t, to)}) ■ (10) 

It is important to mention, that although the delay term 
is skipped in the defining equation for Ao, the delay pa- 
rameters a and t enter indirectly via the trajectory x(t), 
which is inserted in the linearization Eq. ^ . As we show 
later, this dependence is non-negligible and it is possible 
to change the sign of Ao only by variation of a delay pa- 
rameter. 

The consequence of the sign of Ao is documented in 
recent work [T7]. If Ao > 0, we call the resulting chaotic 
dynamics 'strong chaos', else if Ao < but the delay 
system is still chaotic, we define 'weak chaos'. The sign 
of Ao has two major implications. Firstly, for a single 
delay system it determines the scaling of the maximal 
LE with delay time. Secondly, for a network of delay- 
coupled units it determines the possibility to synchronize 
the units at large delay times. If the units are strongly 
chaotic, the delayed coupling cannot compensate for the 
exponential divergence of two nearby trajectories of any 
two systems in the network, and therefore synchroniza- 
tion is not possible. On the other hand, for weak chaos, 
synchronization is in principle possible and depends on 
the network topology. 

In analogy to the general Abarbanel test, we can mea- 
sure the sign of the sub-LE Ao using an auxiliary system 
approach as shown in Fig. [lja) . The system B = TV is a 
copy of the undelayed system A = 1Z. In the presented 



coupling topology both units are receiving the same sig- 
nal A T = T>. If we analyze the equations of motion for 
this setup in the vicinity of the complete synchronization 
manifold (SM) of A and B, we obtain in linear approxi- 
mation the same system, by which we also have defined 
Ao- This means, that in the test setup both units syn- 
chronize if and only if they are in weak chaos. Synchro- 
nization is unstable if and only if they are in strong chaos. 
The result of a numerical simulation of this setup for two 
coupled semiconductor lasers is shown in Fig. [TJb) show- 
ing the cross-correlation between A and B in dependence 
on the coupling strength a. A cross-correlation C = 1 
corresponds to complete synchronization and implies the 
presence of weak chaos. The drop-down of the correla- 
tion indicates transitions from weak to strong and back 
to weak chaos by increasing a. We discuss this behavior 
in detail in Sec. IIVI 



B. 



Scaling of the maximal Lyapunov exponent with 
the delay time 



As we show in the following, the sub-LE Ao, as intro- 
duced by Eqs. (J8j) , determines the scaling behavior of 
the maximal LE A m , as introduced by Eqs. ([3|-([5]), with 
the delay time. Note that there is always a spectrum of 
LEs, in this countable infinite set. A m is the max- 

imal exponent of this spectrum. Further we assume the 
existence of a single attractor and the ergodicity of the 
trajectory x(t) so that we can also skip the dependence 
on the initial conditions. 

We now distinguish between two cases: a) Strong 
chaos, Ao > 0. b) Weak chaos, Ao < 0. In order to 
derive a scaling relation A m (r), we have to assume, that 
Ao(t) = const. In a strict sense this is never true, but 
for sufficiently large r one always observes a saturation 
effect and a decrease of remaining fluctuations in Ao(r), 
so that the assumption of a constant value is valid. 

a. Strong chaos, Ao > 0. We start from the lin- 
earization Eq. Q. If A > 0, for large r the instanta- 
neous term becomes dominant and the delay term be- 
comes negligible. This can be seen from the coordinate 
transformation 



Sx(t) 



which results in 



Sz(t), 



(11) 



Sz= [D J F 1 (x)-A -l](5z + crc" AoT Dr?(x r )^z T . (12) 

For large delay times the delay term in this equation 
becomes exponentially small and can be neglected. The 
resulting Lyapunov exponent of the transformed system 
can therefore be estimated as A z » 0. This implies that 
in the original coordinates A m Ao, meaning that the 
Lyapunov becomes independent of r for large delays. 

b. Weak chaos, Ao < 0. In this case we can es- 
timate the scaling of the LE by considering a strobo- 
scopic sequence of the evolution of our linear system, 
Sx n (0) = 6x(6 + nr) with n E Z and e] - r,0]. One 
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can introduce a growth factor (Lyapunov multiplier) by 



<-n+l 



Mr, 



e)x n ||. so that the LE becomes 



lim - — V In ; 
1-hx It ^— ' 



InC, 



(13) 



where C is the geometric mean of all multipliers. In the 
following we show, that for sufficiently large delay times 
the multipliers /i„ do not depend on t, hence the LE is 
of order r _1 , meaning that in weak chaos a perturba- 
tion of the chaotic system grows on the timescale of the 
delay time. To this end, we introduce the variation-of- 
constants formula, which provides an integral version of 
our initial delay differential equation Q , and evaluate it 
at = 



<5x n+1 (0) = £4(0,-T) 5x n (0) 



a \ dtU n (0,t)DH[x n (t)}dx n (t). 



(14) 



It contains the evolution operator U n (t2, t{) of the auxil- 
iary system Eq. ^ on the n-th r-interval. This operator 
is exponential in t 2 — t\, i.e. with respect to a suitable 
matrix norm it holds that for t 2 > t\ 



||tf(ta,*i)|| <U exp[a(t 2 -h)] 



(15) 



with Ao < a < 0. The bound provided by Uq and a 
should cover potential bursts typical for the linearization 
of a chaotic flow. Because of this exponential bound, 
there is some time tq cx a , such that the term with 
C7(0, — t) can be neglected in Eq. (14), if r > tq. Addi- 



tionally, the integral has only significant contributions in 
a small range close to the end of the integration interval, 
namely for t £ [— tq,0]. This means that a further in- 
crease of r beyond To does not affect the integral. There- 
fore the multiplier introduced above can be estimated by 



|5x n (0)| 



dtU n (0,t) BH [x„(t)] 5x n (i) 



"TO 



(16) 



In leading order, this expression does not depend on the 
delay time. It depends on Ao, a and a set q of other (yet 
unknown) statistical properties of the chaotic trajectory. 
Hence we can write 



A m = -lnC(A ,CJ,q). 
r 



(17) 



For comparison: Given a Floquet problem, in which the 
driving trajectory is r-periodic (x(f) = x(£ + r)) and the 
coupling is linear and diagonal (DH[x(t)] = 1), we obtain 
C = -a/X . 
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Figure 2. (Color online) (a) Maximal Lyapunov exponent 
(LE) A m (solid line) and sub-LE Ao (dashed line) of a single 
laser with self-feedback for a delay time r = 10 ns vs cou- 
pling strength a. (b) Enlarged view of left diagram for small 
coupling strengths cr. 



IV. SINGLE LASER WITH TIME-DELAYED 
SELF-FEEDBACK 

A. Transitions between strong and weak chaos 

We now discuss the most simple system which can show 
strong and weak chaos: a single laser with time-delayed 
self-feedback. The laser dynamics is modeled by the LK 
equations. The parameter values can be found in the 
appendix. If not mentioned differently, all diagrams are 
made for a delay time of t = 10 ns and a pump current 
ofp= 1.02. 

Fig. [2]ja) shows the maximal LE A m (solid line) and 
sub-LE A (dashed line) of a single laser with self- 
feedback in dependence on the coupling strength cr. Al- 
though the coupling strength does not enter directly into 
the conditional equation of the sub-LE, as it would be 
the case e.g. for a linear damping term with a, the sign 
of the sub-LE changes from positive to negative for large 
coupling strengths only through the different dynamics 
x(t) which is inserted into Eq. ([7]). Thus, there is a tran- 
sition from strong to weak chaos with growing a. In the 
region of small coupling strenghts, Fig. |2jb), however, 
there is an additional region of weak chaos where the 
laser is already chaotic, i. e. A m > 0, but the sub-LE 
is still negative. Hence, there is an additional transition 
from weak to strong chaos upon increasing a for small 
coupling strengths. 

Overall, we observe a transition from periodic behavior 
(Goldstone mode with A m = 0) to weak chaos and from 
there transitions to strong chaos and back to weak chaos 
upon monotonically increasing the coupling strength. 

Fig. [3] shows the different behaviors of the maximal 
LE and the sub-LE for the regimes of weak, Fig. |3ja), 
and strong chaos, Fig. [3jb), in dependence on the delay 
time t. As soon as the delay time is large compared 
with the internal time scales of the laser, the sub-LE Ao 
remains constant with increasing r for both weak and 
strong chaos. For weak chaos, A m decreases like 1/r for 
growing r. For strong chaos, A m converges exponentially 
to the positive sub-LE Ao- 

The product A m r is the relevant dimensionless quan- 
tity describing chaoticity in systems with delay. For 
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Figure 3. (Color online) Maximal LE A m (solid line) and sub- 
LE Ao (dashed line) of a single laser with self-feedback for (a) 
weak chaos (a = 21 ns -1 ) and (b) strong chaos (a = 12 ns - ) 
vs delay time r. Product A m r of a single laser with self- 
feedback for (c) weak chaos {a — 21 ns -1 ) and (d) strong 
chaos (a = 12 ns -1 ) vs delay time r. (e) ln(|A m — Xo\/cr) of a 
single laser with self-feedback vs delay time r for strong chaos 
(a = 12ns" 1 ). 



weak chaos, this product saturates at a constant value 
[Fig. [3](c)] which depends on the coupling strength. This 
dependence is caused by different levels of fluctuations 
for different a. Also the delay time r at which the prod- 
uct A m r has reached the saturation value up to a fixed 
distance S depends on a: The closer one is at the critical 
coupling strength cr crit where Ao = 0, the larger the delay 
time has to be, i. e. saturation happens later for larger t. 

For strong chaos, on the contrary, A m r grows linearly 
with increasing r [Fig. J3jd)] since A m gets constant as 
shown in Fig. [3jL). Fig. confirms that the con- 
vergence A m — > Ao > happens exponentially with r. 
The characteristic exponent of this convergence, how- 
ever, is much larger (i. e. the convergence is slower) than 
one would expect from analytic calculations for a simple 
model with constant coefficients. This effect is caused by 
the fluctuations of DF[x(i)] which act like multiplicative 
noise in Eq. @ and Q |2U] . 

Fig. [4] shows the special limit case between strong and 
weak chaos when the sub-LE Ao = for suffiently large 
delay times r. In Fig. |4ja) one can see that the maxi- 
mal LE still decays with increasing r. However, it does 
so very slowly. Fig. |4jb) shows the consequence for the 
product Am r at the critical point: Neither does it grow 
linearly with r like for strong chaos, nor does it saturate 
at a constant value for finite delay times like for weak 
chaos. The point where the product A m r saturates is 
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Figure 4. (Color online) (a) Maximal LE A m (solid line) and 
sub-LE Ao (dashed line) and (b) product A m r (black solid 
line) of a single laser with self-feedback in comparison with the 
least-squares fit of \pr (red dashed line) for critical coupling 
strength where transition between strong and weak chaos oc- 
curs (<r cr it = 13.4 ns" 1 ) vs delay time t. 
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Figure 5. Product A m r of a single laser with self-feedback vs 
coupling strength a for a delay time r = 100 ns. 

only reached for infinitely large delay times r — > oo. A m r 
rather grows like ypr as is shown by the least-squares fit 
with yfr (red dashed line) in Fig.jijb). It turns out that 
this behavior can be explained as well with the effect of 
the multiplicative fluctuations of DF[x(f)]. These fluctu- 
ations can be shown to appear in simple systems in which 
they are modeled as multiplicative white noise [20] . 

Fig. j5ja) shows that for weak chaos and a very large 
delay time r — > oo, the product A m r undergoes a phase 
transition and diverges in the proximity of the critical 
coupling strengths cr cr it,i (gray line) and cr C rit,2 (black 
line). 

We were able to find the occurence of strong and weak 
chaos and the transitions between them by changing the 
coupling strength a also for the Rossler and Lorenz dy- 
namics. The Stuart-Landau, FitzHugh-Nagumo and con- 
tinuous Ikeda dynamics only show weak chaos. 

B. Scaling of the sub-Lyapunov exponent with the 
laser pump current 

We study the dependence of the sub-LE on the laser 
pump current. Fig. |6ja) shows the sub-LE as a func- 
tion of the feedback strength a for different values of the 
pump current. We find that the curves all follow the same 
pattern described in the previous section. The sub-LE is 
negative for small and for large values of the feedback 
strength. The agreement is even quantitative. We find 
that the sub-LE and coupling strength scale with the 
square root of the effective pump current y/p — 1 above 
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the lasing threshold. Fig. |6|b) shows the maximum of 
the sub-LE Ao. max , the feedback strength ao.max for which 
the sub-LE is maximal, and the critical feedback strength 
oo,crit where Ao — as a function of the effective pump 
current p— 1. The data is shown in a double-logarithmic 
scale, indicating a slope of 1/2. In Fig.[6jc) the sub-LEs 
are rescaled with y/p — 1. We find that this scaling rela- 
tion holds well for small coupling strengths. For larger 
coupling strengths the sub-LEs diverge for different pump 
currents. 

To explain this scaling behavior we introduce the fol- 
lowing rescaling of the laser parameters 



N = 



E 



G 



x 



/ r e 



K 



jN so l y/p~l 

1 a 



8= ^TG N N sol y/p^lt. 

Such a scaling reproduces the scaling behavior that we 
found numerically: The coupling strength K and the 
Lyapunov exponent (which scales inversely with time) 
scale with yjp — 1. The LK equations, Eqs. ([!]), can then 
be rewritten as 



dE _ 1 + 

~ds ~~ 2 
dN _ 7 " 
~dJ ~ f 



l a 



1 



N E + Ke ie E(s — y/p-lr) 
N , 



C\/P~ 1 



(18) 



with c = J G^N SQ \jY . We assume small coupling 
strength K <C 1, small carrier densities N = 0(K), and 
a reasonably high pump current c ^/p — 1 = O(l). The 
ratio of photon and carrier life times, 7/T, is a small pa- 
rameter as well. In leading order we obtain the following 
equations: 



dE 

ds 
dN 

ds 



1 



Id 



N E + Ke ie E(s - c yjp - It) 
|£| 2 )- 



(19) 



These equations only depend on the pump current in the 
value of the time delay. Since the exact value of the time 
delay does not influence the sub-LE, we recover the scal- 
ing behavior found numerically. For increasing coupling 
strength K the mapping becomes less exact as can be 
seen in Fig. [6jc) . This results from the fact that the 
rescaled model Eq. ( 19 1 is a weak coupling approxima- 
tion. 



C. Auto-correlations and space-time patterns for 
strong and weak chaos 

In this section we discuss the question whether the dif- 
ference between strong and weak chaos can be identified 
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Figure 6. (Color online) (a) Sub-LE Ao of a single laser with 
self- feedback for a delay time r = 10 ns and several different 
pump currents p vs coupling strength a. (b) Maxima Ao, m ax 
of the sub-LEs Ao, coupling strengths oo.max of these maxima, 
and critical coupling strengths <7 cr it where Ao = for several 
different pump currents p vs effective pump current p — 1 in 
double- logarithmic scale, (c) Data collapse of rescaled sub- 
LEs Ao/Vp — 1 vs rescaled coupling strengths cr/y/p — 1 for 
eleven different pump currents p ranging from p — 1.02 to 
p=1.50. 



from the laser time series itself. Fig. [7] shows two exam- 
ple trajectories (intensity traces) of a single laser with 
self-feedback for the regimes of weak and strong chaos. 
For both strong and weak chaos, there is a characteristic 
structure of high peaks which is significantly higher for 
weak chaos than for strong chaos. This may be caused 
by the larger coupling strength. 

Fig. [8] shows the laser trajectories for weak and strong 
chaos in space-time diagrams where the vertical axis de- 
notes the number of the current delay window of length 
r and the horizontal axis denotes the time offset t in this 
delay window. In such a representation, chaos evolves 
horizontally in space direction for strong chaos. This is 
due to the divergence between two nearby trajectories on 
the internal time scale of the laser which is much shorter 
than the delay time r. In contrast, weak chaos evolves 
vertically in the (discrete) time direction since the diver- 
gence between two nearby trajectories happens on the 
long time scale of the delay. Accordingly, the islands 
of high intensity (red) extend vertically farther in time 
direction for weak chaos than they do for strong chaos. 
In both kinds of visualization (intensity trace and space- 
time diagram), however, one cannot strictly distinguish 
between strong and weak chaos in a nonambiguous qual- 



7 





120 






100 




4 


80 






60 






40 






20 



i;/ 








(a) 



10 

t(ns) 



15 



20 




Figure 7. Example trajectories (intensity traces) of a single 
laser with a self-feedback of r = 10 ns for (a) weak chaos 
(a — 21ns -1 ) and (b) strong chaos (a = 6 ns" 1 ) vs time t. 
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Figure 8. (Color online) Space-time diagram of a laser tra- 
jectory (reddish: high intensity, blueish: low intensity) for 
a single laser with a self-feedback of r — 10 ns for (a) weak 
chaos (a = 21 ns -1 ) and (b) strong chaos (a = 6 ns -1 ) vs the 
number of the current delay window of length r (vertical axis) 
and the time offset t in this window (horizontal axis). 



itative way. 

Fig. [9] depicts the time-shifted auto-correlations C au t 
of a single laser with self-feedback for the regimes of weak 
and strong chaos. For weak chaos, one clearly sees the 
high auto-correlation peaks at multiples of the delay time 
r. Although for strong chaos the chaotic behavior evolves 
predominantly on the internal time scale of the laser, 
whereas for weak chaos it evolves on the time scale of 
the delay, there are non-negligible auto-correlations after 
multiples of the delay time r even for strong chaos. How- 
ever, they are significantly smaller than for weak chaos 
and decay faster with increasing time shift A. 

Considering the auto-correlations C au to.A=T after one 
delay time r in Fig. 10 ^a) , we observe that they do not 
decay for large delay times with increasing r but re- 
main constant. Surprisingly, this is not only the case 





Figure 9. Time-shifted auto-correlations Cauto of a single laser 
with a self-feedback of r = 10 ns for (a) weak chaos (a = 
21 ns _1 ) and (b) strong chaos (a — 6 ns -1 ) vs time shift A in 
units of the delay time r. 



1.0 
0.8 
£.0.6 
1 0.4 
0.2 
0.0, 



(a) 



20 



40 60 
r (ns) 



80 100 



1.0 

0.8 
< 0.6 
° 0.4 
0.2 
0.0, 







(b) 



10 15 20 25 30 



o-(ns _1 ) 



Figure 10. (a) Auto-correlations C au to with a time shift of 
A — t for a single laser with self-feedback for weak chaos (a = 
21 ns -1 , solid line) and strong chaos (a — 6ns~\ dashed line) 
vs delay time r. (b) Auto-correlations C au to with a time shift 
of A = r for a single laser with a self- feedback of r = 10 ns 
vs coupling strength a. 



for weak chaos but also for strong chaos. In Fig. [I0[b) 
we depict the auto-correlations C au to,A=T with a time 
shift of one delay time r in dependence on the coupling 
strength a. We find no sharp transition at the critical 
coupling strengths where the transitions between strong 
and weak chaos appears. For weak chaos, however, the 
auto-correlations C au to,A=r are significantly higher than 
for strong chaos. 

We conclude that a linear measure like the auto- 
correlation function cannot clearly uncover the difference 
between strong and weak chaos. Instead, we propose the 
usage of nonlinear measures in order to detect a more 
significant difference in the relationship between x(i) and 
x(< — t). The idea is motivated by the fact, that in weak 
chaos for sufficiently large delay times, the system can be 
considered to be in a state of generalized synchronization 
with its own time-delayed feedback, whereas in strong 
chaos the state of the system is formally independent of 
its input. Note that this independence does not imply the 
total absence of linear correlations, as we have demon- 
strated. The nature of the nonlinear functional depen- 
dence x[x r ] in weak chaos is exploited in the correspond- 
ing Abarbanel test as described in Sec. |III A] Hence, the 
detection of strong and weak chaos from time series can 
be reduced to the problem of detecting generalized syn- 
chronization |21[ [22] . We refer to the relevant methods 
used in this context, like the evaluation of nearest neigh- 
bors [23], mutual information or transfer entropy |24| . 
However, we assume severe computational difficulties re- 
garding memory usage and runtime connected with the 
necessary delay embedding. Additionally, it might even 
be impossible to detect generalized synchronization from 
a finite time series. 



D. External cavity modes for strong and weak 
chaos 



We can relate the strongly or weakly chaotic behavior 
of a laser subjected to delayed feedback to the properties 
of the external cavity modes (ECMs) of the laser. These 
ECMs are rotating wave solutions of the LK equations 
of the form E{t) — E e lut and n(t) — n with constant 
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amplitude Eq, frequency uj and carrier density n of the 
laser. The spectrum of ECMs is often represented in the 
(uj, n)-plane. In this plane the ECMs lie on an ellipse. 
ECMs located on the lower half are focus solutions called 
modes. The solutions located on the upper half of the 
ellipse are saddle points also referred to as antimodes. 

Depending on the laser parameters, the ECMs have 
different stability properties. In the chaotic regime the 
ECMs can be seen as the skeleton of the chaotic attrac- 
tor |25[ In the low frequency fluctuations regime, 

which occurs for low pump currents and moderate to 
strong coupling strengths, the intensity slowly increases, 
followed by a sudden power dropout. During the buildup 
process, the trajectory travels along the modes in the di- 
rection of the maximal gain mode, until the trajectory 
is expelled along the unstable manifold of an antimode. 
This causes the power dropout. In the coherence col- 
lapse regime, the dynamics can be described as a chaotic 
itinerancy between modes and antimodes. 

The stability can be calculated by inserting the ECM 
solution into Eq. Q. In the long delay limit, the charac- 
teristic equation of a steady state, such as the ECM solu- 
tions, has two types of solutions, which show a different 
scaling behavior with the delay time [27 . The strongly 
unstable spectrum consists of isolated points, which are 
approximated by the unstable eigenvalues of the Jacobian 
of the LK equations without delayed terms. These eigen- 
values do not scale with the delay. Beside this strongly 
unstable spectrum, the characteristic equation has an in- 
finite number of solutions, forming the pseudocontinuous 
spectrum. The real part of these solutions scales inversely 
with the delay. 

One can thus distinguish between strongly and weakly 
unstable ECMs in an analogous way as we distinguish 
between strong and weak chaos. The local eigenvalues of 
the Jacobian without delay terms play a similar role as 
the sub-LE Ao: Strongly unstable ECMs have unstable 
local eigenvalues, and thus a strongly unstable spectrum. 
The maximal eigenvalue is approximated by these local 
eigenvalues and does not change with the delay. The 
weakly unstable (and the stable) ECMs have stable local 
eigenvalues. The strongly unstable spectrum does not 
exist in this case, these ECMs only have a pseudocon- 
tinuous spectrum. Hence, the real part of the maximal 
eigenvalue scales inversely with the delay, just like the 
maximal LE A m for weak chaos [28J. 

In Fig. [TT] we show the projection of the laser dy- 
namics onto the external cavity modes. In the strongly 
chaotic regime, all the modes involved in the dynamics 



are strongly unstable, as illustrated in Fig. 11 'a). Around 
the transition point between strong and weak chaos, a 
few (two or three) weakly unstable modes are involved in 
the dynamics, as shown in Fig. [TT^b) , while most of the 
modes involved in the dynamics are weakly unstable in 
the weakly chaotic regime, as shown in Fig. [Tl|c). These 
features are independent of the pump current. The anti- 
modes are always strongly unstable, both in the weakly 
and strongly chaotic regime. 
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Figure 11. (Color online) Projection of the laser dynamics 
onto the (w,n)-plane for a delay time r = 10 ns and a pump 
current p = 1.10. Thick blue dots denote weakly unstable 
or stable ECMs, thick red dots represent strongly unstable 
ECMs. Panel (a) shows a strongly chaotic trajectory for a = 
20 ns -1 . The modes and antimodes involved in the dynamics 
are strongly unstable. In panel (b) the dynamics around the 
transition point between strong and weak chaos is shown for 
a — 29.28 ns _1 . The trajectory approximates the transition 
point between weakly and strongly unstable modes on the 
ellipse. Panel (c) shows weakly chaotic laser dynamics for 
a = 40 ns -1 . All the modes involved in the dynamics are 



weakly unstable. 



E. Sub-Lyapunov exponents for two self-feedbacks 
with time-scale separated delays 

Until now we have only considered a single laser with 
self-feedback and have discussed the sub-LE in this con- 
text. Now we generalize our investigation by introduc- 
ing a second self-feedback which has a delay time that is 
much smaller than the first one. The linearized equation 
describing the maximal LE is 

Sx = DF(x) (5x+cr s Dif(x T J (5x Ts +(Ti DiJ(x Tl ) 5x Tl (20) 

where t s is the delay time of the shorter self-feedback in 
comparison with the longer self-feedback with delay time 
t\. Additionally to the sub-LE Ao defined from Eq. ^ 
we introduce another sub-LE Ao. s defined from 



Sx 0s = DF(x) (5x , s + 0s BH (x T J <5x 0jf 



(21) 



Thus, for Ao.s we consider a subsystem which includes 
the shorter self-feedback but omits the longer one in the 
linearized equation. As before, the inserted dynamics 
x(t) is the trajectory of the full system which includes 
both the long and the short self- feedbacks. 

Fig. 12 'a) shows the sub-LEs Ao and Ao, s together with 
the maximal LE A m for the case when t\ is much larger 
than t s such that we have a time scale separation be- 
tween the two delays. We observe that Ao s < A holds 
for strong chaos and Ao, s > Ao for weak chaos. Further- 
more, <7 cr it.o ~ Ccrit.o.s holds for the right transition at 
large a, i.e. both sub-LEs change their signs at approx- 
imately the same coupling strength and can hence both 
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Figure 12. (Color online) Maximal LE A m (solid line) and sub- 
LEs Ao,a and Ao (dashed and dotted lines) for a single laser 
with two self- feedbacks n and r s for (a) n = 10 ns, r s = 0.1 ns 
and a B = 5 ns -1 vs coupling strength a :— a\ and (b) r s = 
0.1ns, a B = 5 ns -1 and o\ = 10ns _1 (strong chaos) vs delay 
time r := t\. 



be used as an indicator for strong or weak chaos there. 
For small a, however, only Ao jS changes its sign indicat- 
ing the transitions between weak and strong chaos there. 
Thus, Ao.s is the new relevant sub-LE which determines 
the occurence of strong or weak chaos in a system with 
two time-scale separated delays. This is confirmed by 
Fig. 12 b) which shows that in the regime of strong chaos, 
A m — > Ao.s holds for increasing t\, while for weak chaos, 
A m ~ 1/71. 



V. NETWORKS OF LASERS WITH 
TIME-DELAYED COUPLINGS 



The stability properties of the synchronized state s(t) 
are well-described by the master stability formalism from 
Pecora and Carroll [29 . For completeness, we summa- 
rize the main ideas introduced in their work. A small 
perturbation is applied to the synchronized state, such 
that x z (i) = s(t) + dx. l (t). The equations of motion are 
then linearized around s(t) and the perturbations obey 
the equations of motion 



N 



£x l = DF(s) Sx l 



3 = 1 



(24) 



This set of equations can be decoupled into the ampli- 



tudes £ (t) of the network eigenmodes g € 



for 



which G ■ g — 7a_, g and k = 1 . . . N. With Sx 1 — 
J2k 9 k ' 1 £ k we obtain N decoupled equations 



£ fc =VF(s)Z k + a lk DH(s T )C 



(25) 



Integration of this equation for the fc-th perturbation 
mode g fc yields the maximal LE A&, which tells us 
about the stability of the mode, meaning that a per- 
turbation 5x(t) = (Sx 1 (i),Sx 2 (t),...,Sx N (t)) in direc- 
tion of g k grows or decays exponentially. By construc- 
tion, there exists at least 71 = 1 with the eigenvector 
g 1 = (1, 1, . . . , 1) T , which corresponds to a perturbation 
8x(t) = 8s(t) within the SM. It determines, whether the 
synchronous dynamics is chaotic (Ai > 0) or not. The 
other modes correspond to linear combinations of differ- 
ences between the laser elements and therefore the nec- 
essary condition to find stable complete synchronization 
is X k < for k = 2 . . . N. Since all |7fc| < 1, the knowl- 
edge of the master stability function A(7) is sufficient to 
predict, whether a network of coupled lasers described by 
the coupling matrix G is able to display complete syn- 
chronization or not. 



A. The master stability formalism 

We investigate complete chaotic synchronization of N 
identical coupled laser elements. The state of each ele- 
ment is described by x 2 with i = 1 ... TV. We start from 
the physically motivated ansatz 



N 

x' =F(x i )+a^G^'H(x^). 



(22) 



The matrix G £ R ivxJV contains the network topology. 
G lJ is the normalized coupling strength with which the 
laser i receives input from laser j. In order to guaran- 
tee the existence of the completely synchronized state 



x 1 ^) = x 2 (t) = ... = x N (t) =: s(t), the row sum 
of G has to be constant for all rows, and we choose 
J2j G % i = 1. Then the synchronized state s(t) reduces 
Eq. p2| to 



s = P(s) + o-H(s T 



(23) 



B. Master stability function for weak chaos 

Here we show the general master stability function for 
the limit of large delays in weak chaos. To this end, we 
make use the scaling behavior A m = C(t _1 ). The initial 



point of our considerations is a generalization of Eq. ( 25 1 
which reads 



y = A(t)y + K B(t)y T , 



(26) 



where A(t) := DF 1 [s(t)] and B{t) := T>H[s(t - r)]. The 
exponent A y («;) provides the master stability function. 
Transformation with z = exp(— A y t)y yields a system 
with exponent A z = 

z = (A(t) - A y • 1) z + k c~ A ^ r B(t) z T . (27) 

Weak chaos has the important property, that by changing 
A(t) -> A'(t) = A(t) + e ■ t with e e R and |e| < A , we 
affect the LE only in the order of t _1 . Therefore, if e 
itself is decreasing with r, the effect on the exponent is 
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Figure 13. (Color online) Product A - r of the master stability 
function A(7) and the delay time r = 10 ns for the SM and 
several different coupling strengths a and pump currents p in 
weak chaos regime vs ln(|7|). 



of a smaller order in t, i.e. 0(t p ) with p < — 1. Then 
for sufficiently large r, removing A y from the first term 
on the RHS of Eq. (127} leads to 



y' = A{t)y' + KC - X - T B(t)y' T , 



(28) 



with the exponent Ay = 0(t p ). Comparison of the 
equivalent Eqs. (26 1 and (28 1 shows, that in the lead- 
ing order of t "^reseating of the coupling strength 
k — > k' — k exp(— A y t) led to an exponent A y ~ 0. From 
the knowledge of the zero-crossing at k' , we directly ob- 
tain the scaling law 



M«0 = - In 4- 



(29) 



Although in general the value of k' is unknown, this log- 
arithmic law allows us to connect the maximal LE of the 
SM with the stability of all transversal modes of a net- 
work. We refer to Eq. ( 25 1 , in which k = ajk- Assume 



the maximal LE A m = Ai of the SM is known. It corre- 
sponds to 7! = 1, so A m = Xi = A(<7 ■ 1) is the point at 
which we can fix the master stability function. Then for 



an arbitrary 7^ we obtain from Eq. ( 29 1 the exponent 
M°"7fe) = A m + In |7fe|. 



(30) 



Here we have made use of the fact, that the exponent 
depends only on the absolute value |«| for large de- 
lays. This scaling relation also holds for steady states 
|14j . Transversal stability of the fc-th mode is given if 
A(c7fc) < 0. This leads to the synchronization criterion 



iTfcl < e ~ 



(31) 



which connects transversal and longitudinal stability in 
a network exhibiting weak chaos. 



small in the master stability equation (25 I for r — > 00. 
Hence, any perturbation mode with eigenvalue 7 in a net- 
work, or the mode for complete synchronization in any ar- 
bitrary network, respectively, is unstable for strong chaos 
if the delay time is large. Consequently, complete syn- 
chronization of networks with identical, strongly chaotic 
units is excluded on principle for r — > 00. 

For weak chaos, A(7) basically exhibits a logarithmic 
dependence on 7. However, we observe two deviations 
that can be considered as effects of the finiteness of the 
delay time r which is present in every experiment or sim- 
ulation. As a first deviation, A(7) does not diverge to 
—00 for 7 — > but has a finite value A(0): the sub-LE 
Ao- It is important to note here that the sub-LE Ao is 
equal to the LE of any perturbation mode with 7 = of 
a network (e. g. the mode of complete synchronization 
of the two outer lasers in a bidirectionally coupled chain 
of three lasers) since then the delay term drops out in 
the master stability equation ( p5| , too. This is also the 
case in the experimental setup for the measurement of 
the sub-LE presented in Fig. [TJa). There, the mode of 
complete synchronization has 7 = 0, and hence synchro- 
nization between the two laser happens if and only if they 
show weakly chaotic behavior. 

As a second deviation for weak chaos, Fig. [13] shows 
for several different coupling strengths and pump cur- 
rents that the exponent v in A ~ ln(|7|") = v ■ ln(|7|) 
is not exactly One as for simple analytic models with 
constant coefficients, but rather shows some deviation to 
exponents smaller than One. We were able to find this 
deviation also for the Lorenz dynamics. For constant 
delay time r and constant laser pump current p, the ex- 
ponent gets closer to One the larger a is, i. e. the weaker 
the chaos becomes. Additionally, if the coupling strength 
a is adjusted such that the sub-LE Ao is the same for two 
different pump currents, then the slope is closer to One 
for a small pump current than for a higher pump current. 
As presented in the previous section, a condition for the 
prediction of synchronization in a network, Eq. (31 1, can 
be derived from the behavior of A (7). Due to the em- 
pirically observed slope smaller than One, this condition 
should be refined to 



I72I < e 



-A m r/u 



(32) 



where 72 is the eigenvalue of the coupling matrix G 
with the second largest modulus. This equation allows 
for a more precise practical prediction of synchronization 
in networks of semiconductor lasers, derived from LK 
equation modeling. 



C. Master stability function for the 
Lang-Kobayashi dynamics 

We now investigate the master stability function A (7) 
for the LK equations dynamics. For strong chaos, A(7) 
is constant since the delay term becomes exponentially 



D. Experimental evidence of strong and weak 
chaos in bidirectionally delay-coupled lasers 

In a system with two bidirectionally delay-coupled 
lasers without self- feedback or multiple delays, identi- 
cal chaos synchronization is not stable due to symmetry 
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Figure 14. Cross-correlation function of two bidirectionally 
delay-coupled lasers in the chaotic regime (r = 63 ns). 
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Figure 15. Fiber-based experimental setup of two mutually 
delay-coupled semiconductor lasers. OC\ t 2 stand for the op- 
tical couplers, PC is the polarization controller, Att is the 
variable attenuator, and PD denotes the photodetectors. 



breaking [30 . However, chaos synchronization can still 
exist in the generalized sense [IS]. Here, we give exper- 
imental and numerical evidence that the implications of 
strong and weak chaos regimes also apply to the case of 
generalized chaos synchronization. 

In the context of delay-coupled elements, correlation 
measures can fail to detect synchronization if the num- 
ber of coupled elements is large [22]. However, the cross- 
correlation function is still a good indicator to identify 
generalized synchronization between two delay-coupled 
lasers. We present an example of such a cross-correlation 
function in Fig. [14] which shows distinct peaks at the 
delay time t and its odd multiples. In this figure, the 
large correlation peak at the delay time indicates that the 
lasers are generally synchronized. We argue that gener- 
alized chaos synchronization is only possible, if the lasers 
are operating in the weak chaos regime. In contrast, a 
low correlation peak at a time-lag r is to be expected, if 
the lasers are operating in the strong chaos regime. 

Our fiber-optics-based experimental arrangement is 



schematically shown in Fig. 15 We use two single-mode, 
fiber-pigtailed, discrete-mode semiconductor lasers, emit- 
ting at 1542 nm. The lasers have been hand-selected 
in order to achieve well-matched parameters. The laser 
temperatures and currents are stabilized to an accuracy 
of 0.01 K and 0.01 mA, respectively. The lasers are bi- 
ased at a current of 1.25 Ith, with I t h = 11.7mA being 
the solitary laser threshold current. As shown in Fig. fl5] 
the coupling path includes two 90/10 optical couplers 
OC\ 2i a polarization controller PC, and an optical at- 
tenuator Att. The maximum mutual coupling obtained 
in this experimental arrangement can be estimated to be 
- 40% of the emitted light. The 10% outputs of OC lj2 
are used for detection. 
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Figure 16. (a) Experimentally measured and (b) numerically 
calculated cross-correlations C cr oss with a time shift of A = r 
for a bidirectionally coupled pair of lasers as a function of 
the coupling strength a. The delay time is r = 10 ns in the 
numerics and r = 63 ns in the experiments, corresponding to 
the long delay limit. The experimental coupling strength is 
normalized to the maximum coupling obtained in the setup, 
which is about 40% of the emitted light. 



The values of the cross-correlation between the inten- 
sities emitted by Lasers 1 and 2 for a time shift of r are 
shown in Fig. 16 a) as a function of the coupling strength 
(j. For strong couplings a > 0.3, a region of large cor- 
relation is found. Decreasing the coupling strength to 
0.1 < a < 0.3 results in a sudden decrease of the corre- 
lation. A second region of large correlation can be iden- 
tified for weak couplings a ~ 0.05. The correlation de- 
creases again for the weakest couplings a < 0.05. The 
numerical results, shown in Fig. 16 b) , agree with the 
experimental results. For large couplings, the two cou- 



pled lasers are highly correlated. A distinct region of low 
correlation can be seen for intermediate couplings, while 
the correlation increases again for the weakest coupling 
strengths. The numerical simulations do not reproduce 
the correlation decrease towards zero coupling since spon- 
taneous emission noise sources are not considered. The 
numerical results of the time-shifted cross-correlations for 



the two bidirectionally coupled lasers shown in Fig. 16 b) 
are similar to the auto-correlations of a single laser with 
self-feedback shown in Fig. 10 b). The parameters used 
in the numerical simulations are listed in the appendix. 

Several dynamical states are being observed for vari- 
ation of the coupling along the correlation curve in 
Fig. 16 a). The lasers operate in continuous-wave with 
noisy fluctuations in the absence of coupling. For 
an increasing coupling strength, the lasers follow a 
quasi-periodic route via undamped relaxation oscillations 
reaching a chaotic state at point A. The delay-induced 
dynamics produces a dramatic increase in the laser opti- 
cal linewidth, which increases from a few MHz to several 
GHz due to the coupling. We present in Fig. [T7] the 
typical shape of the optical spectra of the chaotic laser 
for points denoted as B and C in Fig. 16 a). These two 
spectra are qualitatively similar apart from their different 
width, but correspond to significantly different correla- 
tion and synchronization properties. In order to charac- 
terize the optical spectra, we have measured their width, 
defined as the — 20 dB frequency width. The — 20 dB 
width of the optical spectrum is 9 GHz, 20 GHz, 26 GHz, 
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Figure 17. Optical spectra for points denoted as (a) B and 
(b) C in Fig. [Ma). 



and 35 GHz, for points A, B, C, and D, respectively. 

We observe irregular (chaotic) dynamics in the whole 
range from A to D. However, the correlation plot in- 
dicates qualitative transitions of the synchronization be- 
havior between A and B, and B and C, respectively. 
More precisely, the experimental correlation plot shows 
a large correlation around A, and from C to D 7 indicat- 
ing that the lasers are generally synchronized in these 
regions. The degree of synchroniz ability stems from the 
dynamical regime the lasers are operating in. There- 
fore, we can infer that the lasers operate in a weak chaos 
regime in the two high correlation (synchronization) re- 
gions, i. e. for a > 0.3 and a ~ 0.05. In contrast, a region 
of low correlation measure can be seen for intermediate 
couplings a ~ 0.1. This can be associated with a strong 
chaos dynamical regime. The numerical correlation plot 
shows a similar clear distinction between different corre- 
lation regions, with a window of low correlation around 
a ~ 2.5 ns . This low correlation can be linked to a 
strong chaos dynamical regime. 

Our experimental results on two delay-coupled semi- 
conductor lasers presented here support the sequence 
'weak to strong to weak chaos' with an increasing cou- 
pling strength. Even though the master stability function 
cannot be directly applied to the generally synchronized 
solution, the influence of the dynamical regime on the 
synchronization is clearly substantiated by our experi- 
mental and numerical results. 



E. Networks with several distinct sub-Lyapunov 
exponents 

In a network we can define a sub-LE for each indi- 
vidual unit. If the network is not completely symmetric 
then these sub-LEs may differ from each other, even if the 
units by themselves are identical. For example, in a chain 
of three bidirectionally coupled lasers, the middle laser is 
in a different coupling situation than the outer lasers. It 
receives input from two lasers while the outer ones receive 
input only from one laser. It is important to note here 
that the accumulated coupling strength of the inputs, 
however, is constant for each laser of the chain. The oc- 
curence of strong or weak chaos now depends on both 
sub-LEs present in the network. Fig. [l8^a) shows that 
this dependence is simple: The maximal sub-LE of the 
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Figure 18. (Color online) (a) Maximal LE A m (solid line) 
and sub-LEs Ao of the inner laser (dashed line) and the outer 
lasers (dot-dashed line) of a chain of three lasers (see left) 
for a coupling delay time of r = 10 ns vs coupling strength 
a. (b) Maximal LE A m (solid line) and sub-LEs Ao (dashed 
and dot-dashed lines) of pair of two lasers with distinct pump 
currents (see left) for a coupling delay time of r = 10 ns vs 
coupling strength a. 



network determines whether there is strongly or weakly 
chaotic behavior of the complete network's maximal LE. 
If the largest Ao in the network is negative then there 
is weak chaos, if the largest Ao in the network is posi- 
tive, then the maximal LE A m of the network converges 
to the maximal sub-LE of the network. This means that 
potentially additionally present strongly chaotic units do 
not further increase the value to which the complete net- 
work's maximal LE converges for strong chaos. 

Additionally to the case when the network is not com- 
pletely symmetric, also the individual units in the net- 
work may be nonidentical, e. g. in a pair of two bidi- 
rectionally coupled lasers with different pump currents, 
Fig. 18 b). One laser has a pump current p = 1.02 and 
1.50. It is confirmed that even for this 



the other has p ■■ 
case, the largest sub-LE of the network determines the 
maximal LE of the network and that for strong chaos 
the maximal LE of the network follows the largest sub- 
LE. We further observe that, by making the lasers in the 
network nonidentical, we can even induce more transi- 
tions between strong and weak chaos than the transition 
weak/strong/ weak chaos for a single laser or a network of 
identical lasers as shown before. These additional transi- 
tions are interrupted by intervals of a where the maximal 
LE A m decays to Zero, i. e. windows of periodic behavior 
can be found. 

As a consequence, we conclude that complete synchro- 
nization is excluded on principle for arbitrary networks 
which contain at least one strongly chaotic unit if the de- 
lay time is very large. Cluster synchronization, however, 
is still possible. In this case, the weakly chaotic units 
of the network may form one or more clusters. These 
clusters are driven by the strongly chaotic units of the 
network. 
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Figure 19. (Color online) (a) Sub-LEs Ao for several net- 
work patterns where the lasers receive input from one laser 
(see left) in comparison with Ao of a bidirectional triangle of 
lasers for a coupling delay time r = 10 ns vs coupling strength 
a. (b) Sub-LEs Ao of the inner and outer lasers of a chain of 
three lasers in comparison with Ao of a pair of lasers (see 
left) for a coupling delay time t = 10 ns vs coupling strength 
a. (c) Cross-correlation C between the lasers in several net- 
work patterns for a coupling delay time r = 10 ns vs coupling 
strength a. 



F. Sub-Lyapunov exponents for certain network 
patterns 

In the previous section we have shown that different 
number of inputs from other lasers and their arrange- 
ment can have the effect of distinct sub-LEs. This is the 
case even if the lasers by themselves are identical. In 
this section we address the question if the number of in- 
puts is the only criterion which determines the sub-LE 
for otherwise constant system parameters. We would like 
to emphasize once more that the accumulated coupling 
strength of the inputs is constant for all considered net- 
work patterns. 



Fig. 19 a) shows a comparison of the sub-LEs of sev- 
eral networks where each laser receives input from exactly 
one laser. The SM, represented by a single laser with self- 
feedback, and the bidirectionally coupled pair, for which 
the experimental evidence of strong and weak chaos has 
been provided in Sec. |VD| can be seen as limit cases of 
unidirectional rings containing one and two lasers. They 
are compared with a unidirectional triangle which is a 
unidirectional ring of three lasers. The bidirectional pair 
and the unidirectional triangle have in common that they 
do not have an eigenvalue gap. Hence, they cannot syn- 



chronize isochronically for any coupling strength a as 
shown in Fig. [l9[c) for the example of the bidirectional 
pair. Nevertheless, their sub-LEs are exactly the same as 
the one of the SM represented by a single laser with self- 
feedback. Also the sub-LE of Laser B of the configuration 
presented in Fig. [TJa) is equal to the ones mentioned be- 
fore even in the strong chaos regime where the Laser B 
is not synchronized with Laser A. The reason for this 
is that in all cases the lasers receive coherent input from 
a single source which exhibits the specific statistics of a 
laser trajectory. It does not matter if the statistics comes 
from a synchronized laser or an unsynchronized one, as 
long as it is typical for a single laser. This can also be 
nicely seen in Fig. 19 a) by the sub-LE of a bidirectionally 
coupled triangle. If the coupling strength a is not large 
enough to induce synchronization of this triangle config- 
uration, as seen in Fig. 19 [c), then each laser receives in- 



coherently superimposed input from two unsynchronized 
lasers. Hence, the sub-LE of the bidirectional triangle is 
then different from the sub-LEs of the networks with in- 
put from one laser. As soon as the bidirectional triangle 
synchronizes, however, its sub-LE becomes equal to the 
networks with input from one laser. The reason is that 
the signals from the two other laser are superimposed 
coherently and thus equal the signal from one laser. 



Fig. 19 b) shows the two distinct sub-LEs of the inner 
and outer lasers of a chain of three bidirectionally coupled 
lasers in comparison with the sub-LE of the bidirectional 
pair. If the coupling strength a is large enough for the 
outer lasers to synchronize, as seen in Fig. 19 ^c), then the 
inner laser receives the coherent superposition of the sig- 
nals of the outer lasers. Thus, all three lasers in the chain 
effectively receive input from one unsynchronized laser. 
In consequence, the chain can be reduced to a bidirec- 
tional pair of unsynchronized lasers. Indeed, Fig. 19 b) 



confirms that the sub-LEs of the outer and inner lasers of 
the chain become identical and also equal to the sub-LE 
of the bidirectional pair in this regime. 

Fig. |20|a) shows a comparison of the sub-LEs of sev- 
eral network configurations (depicted on the left side) 
in which the lasers receive input from two other lasers. 
Fig. |20jb) shows a comparison of the sub-LEs of several 
networks (depicted on the left side) where the lasers re- 
ceive input from three other lasers. In both diagrams it 
turns out that for small coupling strengths in a partial 
range of the strong chaos regime, the sub-LE indeed only 
depends on the number of inputs and is equal for net- 
works which have the same number of inputs. For larger 
coupling strengths, however, the sub-LEs differ from each 
other. The reason for this is that the two or three inputs 
correspond to an incoherent superposition of two or three 
laser signals. The level of coherence or incoherence de- 
pends on the topology of the subnetwork from which the 
signals originate. These different levels of incoherence 
play a larger role for the laser receiving the signals for 
increasing coupling strength. 

Fig. [20[c) shows a comparison of the sub-LEs of a laser 
which receives unidirectional input from one, two or three 
mutually uncoupled lasers with self-feedback. The infiu- 
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Figure 20. (Color online) (a) Sub-LEs Ao for several network 
patterns where the lasers receive input from two lasers (see 
left) for a coupling delay time t = 10 ns vs coupling strength 
a. (b) Sub-LEs Ao for several network patterns where the 
lasers receive input from three lasers (see left) for a coupling 
delay time r = 10 ns vs coupling strength a. (c) Sub-LEs Ao 
of a laser driven by one, two or three independent laser with 
self-feedback for a coupling delay time r — 10 ns vs coupling 
strength a. 



ence of the topology of the subnetwork from which the 
laser receives the two or three signals is eliminated by the 
fact that the two or three driving lasers are completely 
independent from each other. Also, the driven laser does 
not feed back any signal to the driving lasers. Hence, 
Fig. [20|c) shows the dependence of the sub-LE on the 
number of the completely incoherently superimposed in- 
puts. For increasing number of inputs, the statistics of 
the received summed signal gets less similar to the one of 
a laser and increasingly similar to noise. We observe that 
with increasing number of inputs, the sub-LE gets larger 
for large coupling strengths a and smaller for small a. 



VI. SUMMARY 

In this paper we have extended the investigations on 
strong and weak chaos of [17] by focusing on the dy- 
namics of semiconductor lasers. After describing the LK 
equations and their linearization as a model for the nu- 
merical simulation of a semiconductor laser with time- 
delayed couplings, we have introduced the artificial sub- 
LE Ao and have explained how to determine its sign by 
experiments. 

Strong and weak chaos show different scaling proper- 



ties of the maximal LE with the delay time. The sign 
of the sub-LE Ao distinguishes between strong and weak 
chaos. The transition sequence 'weak to strong chaos and 
back to weak chaos' upon monotonically increasing the 
coupling strength a of a single laser's self- feedback has 
been shown for numerical calculations of the LK equa- 
tions. At the transition between strong and weak chaos, 
the sub-LE vanishes, Ao = 0. At this transition we found 
A m t ~ y/r. Transitions between strong and weak chaos 
by changing a could also be found for the Rossler and 
Lorenz dynamics. 

Counterintuitively, the difference between strong and 
weak chaos is not directly visible from the trajectory al- 
though the difference of the trajectories induces the tran- 
sitions between the two types of chaos. In addition, a 
linear measure like the auto-correlation function cannot 
unambiguously reveal the difference between strong and 
weak chaos, either. Although the auto-correlations af- 
ter one delay time are significantly higher for weak chaos 
than for strong chaos, it was not possible to detect a qual- 
itative difference. But we could relate the trajectories of 
strong and weak chaos to the properties of the external 
cavity modes of the laser. If two time-scale separated 
self-feedbacks are present, the shorter feedback has to be 
taken into account for the definition of a new sub-LE \§ s 
which in this case determines the occurence of strong or 
weak chaos. We have shown that the sub-LE scales with 
the square root of the effective pump current y/p — 1, as 
well in its magnitude as in the position of the critical 
coupling strengths. 

For networks of delay-coupled lasers, we explained us- 
ing the master stability formalism the condition ('yj-l < 
C -A m r £ Qr gtabig chaos synchronization. Hence, synchro- 
nization of any network depends only on the properties 
of a single laser and the eigenvalue gap of the coupling 
matrix. We refined the master stability function for the 
LK dynamics to allow for precise practical prediction 
of synchronization. We provided the first experimen- 
tal evidence of strong and weak chaos in bidirectionally 
delay-coupled lasers which supports the sequence 'weak 
to strong to weak chaos'. For networks with several dis- 
tinct sub-LEs it has been shown that the maximal sub-LE 
of the network determines whether the network's maxi- 
mal LE scales strongly or weakly with increasing delay 
time. As a consequence, complete synchronization of a 
network is excluded for arbitrary networks which contain 
at least one strongly chaotic laser. Finally, we showed 
that the sub-LE of a driven laser depends on the number 
of the incoherently superimposed inputs from desynchro- 
nized input lasers. 
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Appendix A: Parameters for the simulation of the 
Lang-Kobayashi equations 

If not mentioned differently in the text, the constants 
listed in Tab. [Q were used in the simulation of the LK 
equations. 



Appendix B: Critical coupling strengths in 
dependence on the pump current 

Tab. [II] lists the critical coupling strengths for which 
Ao = and at which the transitions from weak to strong 
chaos and from strong to weak chaos appear, in depen- 
dence on the pump current p. 



Table I. Used constants in the simulation of the LK equations. 
Values are taken from |31| . 



Parameter 


Symbol Value 


Linewidth enhancement factor 


a 


5 


Differential optical gain 


Gn 


2.142 x 10 4 s _1 


Laser frequency 


oj 


2tt c/(635nm) 


Pump current relative to Jth 


P 


1.02 


Threshold pump current 






of solitary laser 


Jth 




Carrier decay rate 


7 


0.909 x 10 9 s _1 


Carrier number of solitary laser 


iVsol 


1.707 x 10 s 


Cavity decay rate 


r 


0.357 x lO^s" 1 



Table II. Critical coupling strengths where Ao = and at 
which the transitions from weak to strong chaos ((J cr i t ,i) and 
from strong to weak chaos (o- C rit,2) happen in dependence on 
the pump current p. 



p 


CTcrit,l (nS 1 ) 


CTcrit,2 (nS X ) 


1.02 


0.80 


13.44 


1.05 


0.96 


20.80 


1.10 


1.44 


29.28 


1.15 


1.76 


35.68 


1.20 


2.08 


40.64 


1.25 


2.40 


45.12 


1.30 


2.72 


48.96 


1.35 


2.88 


52.64 


1.40 


3.20 


56.16 


1.45 


3.36 


58.88 


1.50 


3.68 


61.92 
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